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ABSTRACT 

Knowledge of the large-scale vertical velocity is 
required to predict the evolution of the atmospheric 
planetary boundary layer (APBL). Since naval operations are 
often conducted in data sparse regions, Single-station 
assessments of the vertical velocity are necessary. Data to 
evaluate such assessments were obtained from rawinsondes 
taken at San Nicolas Island, California. Vertical velocity 
estimates obtained by vertical integration of the moisture 
budget equation (Q-method) and by the adiabatic method were 
used to initialize an APBL 24-hour prediction model. RMS 
error statistics on predicted inversion height, potential 
temperature, and specific humidity were computed for the 
forecasts and compared to the RMS errors of a persistence 
forecast. Calculation of the vertical velocity by the 
Q-method showed the most promise. However, no 
Single-station assessment method improved on the persistence 


forecasts. 
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Cne aspect of naval operations is battle 


tions ina fixed geographical area. This 


group opera- 


is currently 


referred to as MODLOC. Several scenarios lend themselves to 


MODLCC operations: the Viet-Nam era carrier 
operations, pre-positioning for a projection 


simply a show of force. Training operations, 


based strike 


of power, 


as well, 


OT 


are 


often conducted in fixed areas to avoid hazards to commer- 


Cial shipping and air traffic with weapons firings, 


privacy in operations. 


A second aspect is the need for self sufficiency 


deplcyed force. Any autonomous ability to 


mental prediction within the battle group 


requirements for shorebased support, decreases 


time for perishable information and reduces the 


communication assets. 


A third aspect is the effect of the atmosphere 


perfcrmance of weapons and sensors that depend 


magnetic and electro-optical wave propagation. 


cal factors which affect such prepagation 


Significantly in four to six hours. The importance of APBL 


prediction to naval operations lies in: 


OL £0 


of a 


make environ- 


reduces 


on 


the 


delivery 


demand on 


+he 


on electro- 


Meteorologi- 


can change 





e Prediction of the existence and evolution of electro- 
magnetic ducts. : 


e Prediction of the inversion height, _ where it is_ known 
that optical propagation is degraded due to turbulence. 


This thesis is an examination of the accuracy of single- 
station assessments of the large-scale vertical velocity, 
and the effectiveness of these assessments for an existing 
APBL prediction model. Three methods of computing the ver- 
tical velocity will be discussed. Vertical velocities cal- 
culated by single-station methods will be used in an APBL 


model to ascertain the relative merits of these methods. 


10 





II. BACKGROUND 


The single-station forecasting problem probably has been 
of interest as a military requirement since the origin of 
armed conflict. Oliver and Oliver (1945) put single-station 
assessments into a naval perspective: 


During the last few years, wartime conditions have made it 
necessary for isclated combat units to issue forecasts in 
regions where no network of meteorological stations could 
be available. Frequently the data from several stations 
or from reconnaissance planes are available, but in some 
regions. the forecaster must rely only on, surface. and 
upper-air observations made at his own Station. This is 
particularly true in the case cf ships at sea. Hence, it 
1s important to Setee ce proficiency at oa Eeeeend informa- 
ticn with limited aerological data.... He must further be 
able to distinguish from solutions which are internally 
consistent, and those which are inconsistent and therefore 
to be discarded. 


In APBL assessments the assumption of horizontal homo- 
geneity (one-dimensionality) is often used. All APBL pre- 
diction models require, as an externally known parameter, 
the large-scale vertical velocity which direcly affects 
changes in the height of the APBL. Furthermore, the height 
of the APBL affects other properties of the boundary layer 
such as the temperature and humidity. 

Tata for this examination were obtained from rawinsondes 
taken at San Nicolas Island (SNI), California. SNI data are 


sufficiently representative of a marine environment for this 
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purpcse. Battalino et al. (1979) have stated that condi- 
tions at SNI may be marine, continental, or mixed depending 
on the circulaticn patterns. A decision to use the SNI data 
was made, in part, because of the availability of siailar 
data from Vandenberg (VND) Air Force Base, California and 
Point Mugu (PMTC), California, so that examination of the 
horizontal variability of the area and divergence derived 
vertical velocity values could be accomplished. Fags 1 is 
provided for geographic orientation. 


1a Ww 11h W 11 W 117° Ww 
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Figure 1. Southern California Ccastal Waters. (After 
Battalino et al. 1979.) 
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A. THE PREDICTION MODEL 

An APBL prediction model in which the mean vertical 
moticn is an external parameter is used to establish the 
affect of the derived values of vertical “velocity on boun- 
dary layer evolution. The model used is based on the one 
presented by Stage and Businger (1981) and was coded for 


operational use by Brower (1982) . The predictive equations 


ares 
dq/fdt = h-t(wg'(0) + Wag) (1a) 
dQ, /dt = h-1(w'@! (0) + RQ - R, + We. Ad) (1b) 
dh/dt = Wy + We (1c) 


Where q is the total water content (liquid and vapor), h is 


the inversion height, 8 is the equivalent potential temper- 





ature, w'q'(0Q) is the surface flux of moisture, w'@! (0) is 
the surface buoyancy flux, A is the difference between the 
variable value immediately above the APBL and its well mixed 
value in the APBL, R. 1s the rate of radiative heat gain per 
unit area at the cloud base, and R, is the rate of radiative 
heat loss per unit area near the cloud top. We is the 
entrainment velocity and W, is the large-scale vertical 


velocity at the inversion height. 
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The model atmosphere is two layered, composed of a tur- 
bulent lower layer in which equivalent potential tempera- 
ture, specific humidity and wind speed are constant with 
height, and a stable upper atmosphere with constant specific 
humidity and buoyancy gradients. The layers are separated 
at the inversion height by a zero-order discontinuity 
(jump). Fig. 2 depicts how a sample sounding is character- 
ized in the model atmosphere. 

Surface fluxes of moisture and bouyancy are diagnosed 
from bulk formulae (Lilly 1968). The liquid water profile 
is ccmputed adiabatically and reduced by about 30 percent to 
agree with empirical results (Fairall et al. 1981.) Long 
wave radiative cooling is treated using the Stephan-Boltzman 
law, and short wave radiative warming is treated using the 
delta-Eddington approximation (Pairall et al. 1981). 

The large-scale vertical velocity, W, , is treated as an 
external parameter, and is constant over the model run 
(24-hr). The system is closed with the specification of the 
entrainment velocity which is accomplished by the assumption 
that dissipation of turbulent kinetic energy (TKE) is a 
fixed fraction of the production of TKE (Stage and Businger 


1981). A 30-minute time step is used. 
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Tit. METHODS 


At least three methods which can be used to compute the 

large-scale vertical velocity are: 

e The kinematic method. 

e The adiakatic method. 

e Integration of the moisture budget equation. 
The kinematic method requires wind sounding data from at 
least three appropriately spaced stations. The adiabatic 
methcd and the integration of the moisture budget equation 


are truly single-station assessments. 


A. THE KINEMATIC METHOD 

The kinematic method can be used to compute the vertical 
velocity when nearly simultaneous (within two hours) sound- 
ings are available at three or more appropriately spaced 
stations such as VND, PMTC, and SNI (Fig. 1). Following 
Saucier (1955), the vertical velocity can be derived: 

Ow/OzZz = - (Qu/dx + av/oy), (2) 

where u, v and #w are the zonal, meridional and vertical wind 
compcnents. The horizontal wind components can be linearly 


apprcximated by a Taylor series expansion as: 


ue + (QU/dX) X + (QU/dY) Y, (3a) 


or 
ll 


v=zM + (AV/0X) x + (dV/dyY) Y- (3b) 
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Fq. (3) suggests that the horizontal derivatives may be 
obtained from a regression of the horizontal wind components 
onto the horizontal plane (Fairall et al. 1981). Decomposi- 
tion of the wind vector from the sounding data at the three 
stations (SNI, VND, PMTC) into zonal and meridional compo- 
nents is required. Vertical profiles of the horizontal wind 
compecnents were constructed at a 200 meter vertical incre- 
ment by linear interpolation from the irregulariy spaced 
data levels tc those levels that are an integer multiple of 
200 meters. Multiple linear regression of each component of 
the horizontal wind onto the horizontal plane is performed 


as fcllows: 


u at bx + cy, (4a) 


v= 4+ ex + fy, (4b) 

where a through f are the regression coefficients. 
Comparison of Eq. (3) with Eq. (4) indicates that the 

coefficients b and £f give the zonal derivative of zonal wind 

and the meridional derivative of meridional wind at the cen- 

ter cf a triangle formed by the three locations and at each 

vertical level where the regression is done. 


Vertical integration of Eq. (2) is performed to give a 


vertical profile of vertical velocity to three kilometres: 


3km 


w(z) = - So((b) + (f) ) az, (5) 
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where w(z=0) = 0 and Az = 200m. 


Be THE ADIAPATIC METHOD 

For isentropic motions there is no time change of poten- 
tial temperature following the motion. Based on this 
assumption Saucier (1955) specifies the adiabatic method for 
estimating vertical velocity as: 


(d9/dt) + V- (V8) 
ee (6) 
(08/3Zz) 


where 6 is the potential temperature. Saucier further 
states that the local temperature change due to advection 
may ke on the order of that due to vertical motion. In this 
study, calculations of vertical velocity are made both 
including and neglecting the advective term. Description of 
the advection calculation is left to a subsequent section. 
The vertical gradient is determined by a linear fit to 
the lapse rate of potential temperature above the inversion. 
The local time change of temperature is obtained at fixed 
heights (z = 1.0, 1.25, 1.5 km) from sequential soundings at 
the same location. These heights are above the inversion. 
It is assumed that the vertical velocity at the surface is 
zero. The vertical velocity at the inversion is taken as 


the average of the vertical velocity at the fixed heights 
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linearly scaled to the inversion height. Thus, a single- 


station assessment of the vertical velocity is obtained. 


C. INTEGRATION OF THE MOISTURE BUDGET 

Integration of the moisture budget equation to compute 
the vertical velocity was performed by Lenschow (1973). 
This method will hereafter be referred to as the Q-method, 
and is based on the assumption of well-mixed specific humid- 
ity in the boundary layer, and that changes occur due to 
fluxes at the sea surface and inversion only. 

Vertical integration cf the moisture budget equation 
through the boundary layer and across the interface between 
boundary layer andthe stable atmosphere above yields the 
follcwing expression for vertical velocity, w, at the level 
of the inversion, h: 


ce A(qtl) (dh7dt) - h(d(q+1)7dt) + w'q'(0) - 
a(qtl) + (Bh/2) : 





where q and 1 are vapor and liquid water contents, w'q' is 
the moisture flux, 8B is the vertical gradient of (q+l) in 
the mixed layer, and ,A(qtl) is the difference between total 
moisture in the mixed layer and total moisture immediately 
above the mixed layer. With no liquid water, z<<h, and 


neglecting moisture advection Eq. (7) is simplified to: 
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Aq(dh/dt) - h(dgq/dt) + w'q' (0) 
w= ——__ = . (8) 


Lilly (1968) assumes that 8 is zero. Lenschow (1973) found 
Bh/2 to be small and certainly much less than the moisture 
jump. 

Eq. (8) provides another independent assessment of the 
vertical velocity at the inversion. The time derivatives 
(dh/dt and dgq/dt) are determined from sequential soundings. 
The moisture jump and inversion height are taken as averages 
for the two soundings, and the surface moisture flux is det- 


ermined by the bulk method, i-e., 


w'q'(0) = Ch Vig (F5 a Tio) (9) 
The exchange coefficient, C,, is computed from the wind 


speed dependent relation of Large and Pond (1980). The sur- 
face specific humidity, q,, is determined via Tetens' for- 


mula for the vapor pressure (Buck 1981). 


De DETERMINATION OF ADVECTION FROM THE THERMAL WIND 

Under the geostrophic assumption, the horizontal temper- 
ature gradient is directly proportional to the vertical 
shear of geostrophic wind. Following Byers (1974): 


v= £-1(3¢/8x), 3 uw = -£71(3¢6/dy),, (10) 
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where v and u are the meridional and zonal components of the 
geostrophic wind, f is the Coriolis parameter, @? is the 


geopctential and the horizontal derivatives are on the con- 


stant pressure surface. Vertical differentiation of 
Eq. (10) and substitution from the hydrostatic equation 
yields: 

Ov/op = f-!(-a) , (11a) 

du/op = f-*(a) , (11b) 
where a is the specific volume. Substitution from the 


equation of state and rearranging terms, the components of 


the temperature gradient are: 


(OT/dX), -(fp/R ) (dv/dp) -« (12a) 


(OT/dyY), (fp/R )(Qu/dp) « (12b) 
where R is the gas constant for dry ‘air. 

With a wind sounding, assumed to be geostrophic above 
the boundary layer, the horizontal temperature gradient can 
be specified and a single-station assessment of the thermal 
advection *erm in Eq. (6) is obtained: -(¥-VT)=-(V-V@), at 
a constant level. 

This calculation is performed by subtracting the wind at 
the first level above the inversion from the wind at an 


upper level. The advecting wind is the mean wind in the 


layer. 
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Ee SINGLE-STATION ASSESSMENT SCORES 

The separate vertical velocities computed from single- 
station assessments by the adiabatic method and the Q-method 
were both used as the large-scale vertical velocity in the 
APBL model. 

Foth methods required use of finite differences in time 
to approximate time derivatives. The time increment used 
was the time difference between subsequent soundings at SNI; 
it varied from four to eighteen hours. To establish the 
vertical velccity for a model run, the time weighted average 
vertical velocity was used, as computed from the SNI data 
for the 24 hour period following model initialization. 

The model requires initial values of inversion height, 
and initial profiles of potential temperature and specific 
humidity. Sea surface temperature, required for computing 
surface moisture flux, was constant through the model day. 
Wind speed, required for computing surface moisture and 
bouyancy fluxes, was input as a linearly changing function 
of time. 

Three model runs were made on each set of initial condi- 
tions: one with vertical velocity by the Q-method, a second 


with vertical velocity by the adiabatic method and the last 
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with vertical velocity set to zero. The advective term of 
Eq. (6) it was found to be much less than the local time 
change of potential temperature. Thus, it was neglected. 

Root Mean Square (RMS) error was computed for predicted 
inversion height, h, potential temperature, 0, and specific 
humidity, q. 

RMS error = . x = yi ie ee , 

where Nis the number of observations, x, is the observed 
value cf a variable, and y, is the predicted value of a 
variable. Observed values were obtained from rawinsonde 
data available within the 24 hours following the time of the 
initial data. Typically, two or three soundings were avail- 
able for each 24-hour period following initialization. 
Predicted values were taken from the model output for the 


time correspcnding to that of the observed values (+ 15 


minutes, due to the 30 minute model time step). 
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Two periods (cases) were considered. Case I was from 31 
August to 4 September 1981 and consisted of the analyses of 
nine model predictions verified against twenty observations. 
Case II was from 15 to 18 September 1981 and consisted of 
the analyses of eight predictions verified against seventeen 


observations. 


A. CASE I (31 AUG - 4 SEPT 1981) 

Throughout the period the eastern North Pacific Ocean 
was under a high pressure system that was moving slowly 
eastward. Southern California was dominated by a thermal 
trough, leading to weak offshore flow at the surface, while 
California was to the south of the strongest westerlies at 
upper levels. Pigs. 3 through 6 show the National Meteoro- 
logical Center (NMC) surface and 500 millibar (mb) analyses 
for Case I. 

Selected prints of the GOES WEST Infra-red (IR) imagery 
indicate low broken to overcast cloudiness in the vicinity 
SeeoNL (Fig. 7). Surface observations at SNI indicated 


overcast conditions with fog at night, and scattered low 
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cloudiness with haze persisting after the late morning 
burn-off of the fog. ; 

A Pacific high pressure system dominated the eastern 
ocean area west of California on 31 August (Fig. 3). A low 
pressure system was located in the Gulf of Alaska with its 
associated weak frontal systen. Tropical Storm Irwin was 
dissipating southwest of the Baja Peninsula. A thermal 
trough dominated the central California valley. Very weak 
westerly flow aloft was predominant over southern California 
and most of the eastern Pacific Ocean. 

After 31 August the Pacific high intensified, while the 
low pressure system in the Gulf of Alaska moved into the 
northwest United States and the assosiated cold front dissi- 
pated. The thermal trough over inland California moved 
Slightly eastward and dominated the high desert region of 
Arizcna and southeastern California (Figs. 4 through 6). 
Flow aloft remained very weak over southern California and 
the low level flow in the vicinity of SNI was very weak (3 


to 5 kt) and predominantly offshore. 
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30 





B. CASE II (15-18 SEPT 1981) 

Throughout the period the polar front was located in the 
eastern part of the North Pacific Ocean and Gulf of Alaska. 
A high pressure system dominated the area southeast of the 
polar front. A thermal trough was located over the inland 
California valley and through much of the southwestern 
United States and Mexico. Figs. 9 through 13 show the NMC 
surface and 500 mb analyses for Case II. 

Satellite imagery indicates considerable clear area in 
ere) vicinity of SNI. Surface observations at SNI reported 
scattered low and high clouds with occasional reduced visi- 
bility due to fog and haze. Selected prints of the GOES 
WEST IR imagery are shown in Fig. 14 . 

Cn 14 September high pressure prevailed over the eastern 
Pacific Ocean. The thermal trough dominated the southwest- 
ern United States and the Mexican Plateau. A low pressure 
system with a weak frontal system existed in the Gulf of 
Alaska. Weak unorganized flow existed south of the strong- 
est westerlies in a large region of very flat height gradi- 
ents over most of California and the eastern Pacific Ocean 


(Fig. 9). 
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Ey 15 September a large amplitude ridge aloft had built 

c 
over the western United States and Canada (Fig. 10). The 
thermal trough weakened in central California as the ridge 
axis in the eastern Pacific Ocean migrated eastward. Low 
level flow in southern California was offshore and stronger 


(12 to 22 kt) than in Case I. 


C. CIFPERENCES BETWEEN CASES I AND ITI 

Comparison of the mixed layer and overlying data clearly 
indicates that different synoptic situations existed LOT 
the two cases. These differences are summarized in Table I 
and illustrated in Fig. 8 : 


TABLE I 


Comparison of the two cases 


Case [I Case II 
mean h 604 o 262 om 
range of h 511 to 667 on 145 to 357 o 
mean 9 Ga 2: e 13.1 C 
mean A®@ SO 5c 4H4.1C 
mean g 8.3 g/kg 9.8 g/kg 
mean Ag -6.0 g/kg “4.4 g/kg 
mean 09/dz above h 4.2 C/km 228 C/7Kn 
méan 9q/dz above h -.35 g/(kg-km) -.83 g/(kg-kn) 


Cn the average, inversion height was nearly 2.5 times 
greater in Case I than in Case II. Case II had a warmer and 


more moist boundary layer with smaller magnitudes of jumps 
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in beth potential temperature aad specific humidity. The 
vertical gradient of potential temperature (stability) above 
*+he inversion is 1.5 times greater in Case I than Case Il, 
and the vertical moisture gradient is nearly 2.5 times 


larger in Case I than Case If. 


height (km) 





g/kg C 


Figure 8. ge Model Atmosphere: Cas@2 I vs. Case If 


a Case I, dash: Case ITI 


Avera 
Solid 
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Figure 14. GOES WEST IR, Case II. 
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V. RESULTS 


Tirect measurement procedures do not exist to obtain 
vertical velocity values. Therefore an indirect evaluation 
procedure was nesessary. This was accomplished by using the 
RMS Error (section III,E.) for model predictions of inver- 
sion height, potential temperature and specific humidity, 
with the large-scale vertical velocity specified by the 
adiakatic method, Q-method, vertical velocity equal to zero, 
and a persistence forecast. It is noted that the Q-method 
is kased on measured values of the prediction variables 
(Eq. (1)). Therefore the model is expected to perform well 
with peer iSal velocities cemputed in this way. 

In both cases the inversion persisted for several days. 
Eqe (1c) illustrates that W acts to keep a lid on the boun- 
dary layer only when subsidence (W <0) is occurring since W 
provides only for increasing the inversion height, i.e., 
"detrainment™ is not allowed. Then from a pragmatic view, 
negative vertical velocities (subsidence) are sought while 
positive ones are not. Large-scale subsidence is also 
expected, frcem a synoptic viewpoint in both cases, since SNI 


is under the summertime subtropical ridge. 
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Table II gives the averages of the vertical velocities 
for €ach case calculated by the Q-method, the adiabatic 
methcd, and the kinematic method. The standard deviation is 
given in parentheses. Figs. 15 and 16 show the frequency of 
occurrence of calculated vertical velocity for Case I and 
Case II, respectively. Assuming that subsidence existed at 
the inversion height, the adiabatic and kinematic methods 
are unable to diagnose subsidence. The Q-method consis- 
tently gives subsidence where subsidence is expected. Dis- 
cussion of possible sources of errors in these methods is 
left to a subsequent section. 


TABLE IT 


Vertical velocity (cm/sec) comparisons 
Q-method Adiabatic method Kinematic method 


cast 28 GT} 238 WF 3 


A. LISCUSSION OF CASE I 
The RMS errors for Case I appear in Table III. The 
average time from initialization to verification for this 


case was 15.6 hours. Persistence forecasts yield the least 
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RMS error for all variables (h, 9, and q). When the model 
is initialized with the large-scale vertical velocity con- 
puted by the adiabatic method or when the vertical velocity 
is assumed to be zero, very poor results are obtained in the 
inversion height prediction. Vertical velocities computed 
by the Q-method result in better prediction of inversion 
height than the adiabatic method. 


TABLE IIt 
RMS Error (Case [) 


h 9 q 
Persistence 61 o Dee -6 g/kg 
W=0 167 mo Zee € -8 g/kg 
pa ec 106 o 2-4 C «9 g/kg 
Adiabatic method 249 2e2C -7 g/kg 


The large errors in the inversion height for the W=0 
prediction occur because the boundary layer is constrained 
to growth by entrainment. The poor results of the adiabatic 
methcd are also to be expected because of two basic problems 
with this method. Pics it, examination of Eq. (6) reveals 
that the lapse rate is assumed to be constant and the numer- 
ator is approximated by a finite difference in time over six 
to twelve hours. An uncertainty of + one degree in poten- 
tial temperature can lead to an uncertainty in vertical vel- 
ocity of +.5 (cm/sec), , assuming nominal values of the free 


atmosphere lapse rate (2.5 to 5 C/Km). Since the actual 
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vertical velocity may be positive or negative and often on 
the crder .5 (cm/sec), the adiabatic method is not adequate 
to specify either the magnitude or the sign of the vertical 
velocity. Second, one must choose an appropriate time 
interval over which the finite difference applies. In prac- 
tice the time interval was that of the rawinsonde schedule; 
however, the soundings must be scheduled suffiently close in 
time to reveal changes insign of the tendencies of the 
boundary layer variables. 

The Q-method also suffers from the disadvantages of 
finite difference in the specific humidity (dq/dt) and the 
inversion height (dh/dt). However, the moisture flux term 
has no. such problen. In addition, specification of the 
moisture and inversion height variation with time required 
for these calculations does not require profiling devices 
such as radiosondes, but may be obtained from surface 


instrumentation. 


Be CISCUSSICN OF CASE ITI 

Average time from initialization to verification was 
14.6 hours. Table IV summarizes the results. 

None of the predictions did particularly well with 


respect to RMS error of inversion height. The same ordering 
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TABLE IV 
RMS Error (Case ITI) 


h 8 
Persistence 105 leer ale xe 
W=0 375 on fora G 1.49/kg 
_Q-method 142 5 12 2.49/kq 
Adiabatic method 527 Ugicnle 1.3g/kg 


of RMS errors for inversion height prediction occurs with 
+he vertical velocity specification methods in this case. A 
major difference in these results is that the model pre- 
diction of specific humidity with vertical velocity by the 
O-meéethod exhibits a 24% error, much greater than any other. 
This is to be expected since the average inversion height 
was 262m, over 40% lower than in Case I, and in Eq. (1a) the 
specific humidity tendency is inversely proportional to 
inversion height. Hence, low initial inversion heights, as 
was the case, coupled with predictions of inversion heights, 
which are too low, contribute to the over prediction of spe- 
cific humidity. The Q-method yields the least RMS error in 


potential temperature predictions. 


c. CLISCUSSICN OF THE KINEMATIC METHOD 

Vertical velocities calculated by the kinematic method 
were positive more often than not as indicated in Figs. 15 
and 16. They were much larger in magnitude than those cal- 


culated by either the adiabatic method or the Q-method. 
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Therefore application of these to the model would have pro- 
duced larger RMS errors than any presented. Two reasons for 


the poor results are: 


e Improper specification of the bottom boundary condi- 
a0 Nie 


e Small sample size for the regression approach. 

Cver the ocean the bottom boundary condition for Eq. (2) 
is that the vertical velocity equals zero at the sea sur- 
face. This is a good assumption since the surface is flat 
with respect to a local earth tangent plane. However, since 
SNI, PMTC, and VND are coastal stations, the boundary condi- 
tion for Eq. (2) requires specification of the terrain 


(Dutton 1976). 


D. CONC LUSION 

Since both Case I and Case II represent small samples of 
observations the error statistics have little stability. 
However, they do suggest that vertical velocity calculated 
by the Q-method has the most merit as a _ single-station 
assessment. Clearly, the assumption of no vertical motion 
or vertical motion derived from the adiabatic method has 
questionable value in APBL assessments. Fig. 17 shows a 


model prediction with vertical velocity computed by the 
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Q-method. Fig. 18 is the same model run but with upward 
vertical velocity computed by the adiabatic method. Compar- 
ison of these figures shows the good agreement in gq pre- 
diction indicated in Table III . On the other hand, the 
superior result of the Q-method is well illustrated by con- 
parison of the evolution of the inversion height. The 
importance of diagnosing subsidence when it exists is evi- 
dent by this comparison. 

The assumption that SNI represents a marine environment 
is a possible source of error in the Q-method. Also, sea 
surface temperature has an indirect effect on the calcula- 
tion of vertical velocity through Eq. (9). The sea surface 
temperatures used were obtained from the Sea Surface Thermal 
Analyses of the National Oceanographic and Atmospheric 
Administration, rather than being in-situ measurements. 

Additional research is needed in specification of the 
large- scale forcing of the APBL. Vertical velocity by the 
Q-method appears promising, however larger samples must be 
obtained to gain statistical stability. This need not be an 
expensive venture. One only needs continuous measurements 
of sea surface temperature, air temperature, humidity and 


Wind. That can be done on a near-continuous basis with 
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currently available instrumentation that does not place 


great demands on personnel. 
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Pigure 17. Model prediction: vertical velocity by Q-method, 
Initialization: 0900 Pacific Daylight Time, 
1 September 1981. SE Coane well-mixed 
Specific humidity (solid), Jump (dotted). 
Middle frame: well-mixed potential temperature 
(solid), jump (dotted Bottom Frame: 


Inversion height inteticayis lifting condensation 


level (dottedj. bserved data are shown with an 
asterisk (@). 
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Figure 18. Same as 17 except Adiabatic Method 
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